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O ■ ABSTRACT 

O 
(N 

^ I We investigate the effect of X-ray ecfio emission in gamma-ray bursts (GRBs). 

y^ \ We find tliat tfie ecfio emission can provide an alternative way of understanding 

■ \ X-ray sliallow decays and jet breaks. In particular, a sliallow decay followed by 

CSJ \ a "normal" decay and a further rapid decay of X-ray afterglows can be together 

explained as being due to the echo from prompt X-ray emission scattered by dust 

grains in a massive wind bubble around a GRB progenitor. We also introduce an 

extra temporal break in the X-ray echo emission. By fitting the afterglow light 



O 

^ I curves, we can measure the locations of the massive wind bubbles, which will 

j^ \ bring us closer to finding the mass loss rate, wind velocity, and the age of the 

progenitors prior to the GRB explosions. 
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O ' Subject headings: dust, extinction — gamma rays: bursts — ISM: bubbles 
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scattering — stars: mass loss — X-rays: general 



I> ; 1. INTRODUCTION 

o 

Gamma-ray bursts (GRBs) are the most luminous and intriguing explosions in the uni- 
/\^ • verse. A number of physical interpretations for GRBs have been proposed since they were 

c3 ■ discovered about four decades ago. The most suc cessful model for GRBs and their after- 



glows is the fireball shock model (JMeszarod l2006l . and references therein). Although the 



interpretation for the prompt gamma-ray emission is still controversial even in this stan- 
dard scenario, the long-wavelength afterglows have been w idely accepted to be the emi ssion 



from relativistic shocks sweeping up an external medium (JMeszarosi l2006l : IZhangI 120071 . and 
references therein). 
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In spite of the success of the standard external shock model for interpreting most features 
of afterglows, there a re still a few issues recently raised by the Swift X-ray Telescope (XRT) 
( JGehrels et al.l l2004f). The first issue is the ubiqu i tously detected shallow decay of early 
X-ray afterglows (JNousek et al.l l2006l : IZhang et al.l l2006l ). In general, to account for the 
shallow decay, a "ref reshed shock" scenario is assumed, either by a conti nuous activity of 
the GRB progenitor JPai fc LulEoQsl : Izhang. fc MeszarosI I2OO1I: [Pai 20oi ) or by a power- 
law d istribution of the Lorentz factors in the ejecta (JRees fc Meszaroslll998l : ISari &: Meszaros 
20001 ). However, although the "refreshed shock" model is generally able to explain the flatter 
temporal X-ray slopes from current analysis, it may not be able to explain self-consistently 
both the X-ray an d the optical aftergl ows of some well-monitored GRBs, such as GRB 050319 
and GRB 050401 JPan fc Piranlbood ). 



The second issue is the lack of jet breaks that are supposed to be detected simultaneously 
in X- ray and optical light curves (jPanaitescu et al.ll2006l : ISato et al.ll2007l : lBurrows fc Racusin 
20071 ). This raises a serious concern about our understanding of the afterglow emission. Cur- 
rently, possible solutions to this pro blem are additional mechanisms taking place in the early 
X-ray emission flFan fc Piranll2006 ) , X-r ay emission arising from different outflows with the 



optical emission (jPanaitescu et al. 



20061 ). or alternati vely, the tempor al break in X-rays is 



masked by some additional source of X-ray emission (ISato et al.l 120071 ). so that a jet break 



at the expected time could still be seen in the optical band. T 



bring renaissance to the physical meaning of a jet break ( JRhoads 



le la st interpre tation may 



2OOII ). and therefore revive the relation between Ej^t and -Epeak (JGhirlanda et al.ll2004j ). 



1999h and E;^, [^Frail et al. 



Recently IShao and Pail (120071 ) realized that the impulsive prompt X-ray emission scat- 
tering off the surrounding medium could give rise to a long-term X-ray flux, which is compa- 
rable to that of most observed afterglows. The scattered X-ray emission (which is called echo 
emission herein) could b e an additiona l contributor to the early X-r ay afterglow, which was 
or igin ally p ut forward bv iKlose I (119981) and then furthe r discu ssed by iMeszaros fc Gruzinov 



( I2OOOI ) and ISazonov &: Sunyaev I (120031 ). IShao and Pail (120071 ) also worked out detailed light 
curves of the echo emission that could nicely reproduce most of the temporal features of 
observed X-r ay afterglows, i.e., a shallow decay fol lowed by a "normal" decay and a further 
rapid decay (jZhang et al.ll2006l : iNousek et al.ll2006l ). This scenario provides a self-consistent 
solution to both the shallow decay and the jet break proble m. By fitting the X-ray light 
curves of GRB 060813 and GRB 060814. IShao and Pail (120071 ) concluded that the scattering 
medium, i.e., a dusty shell should be at a distance of about tens of parsecs from the GRB 
source. 



Coincidentally, this conjectured dusty shell by IShao and Pail (120071 ) is consistent with 
the massive wind bubble around a Wolf-Rayet (WR) star, which has been proposed as the 
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progenitor of a GRB. The study of stellar structure and evolution has provided a mass- 
loss history of a massive star, starting with the Main Sequence (MS) stage as an O star, 
which will evolve into a WR star and then produce a final supernova and/or a GRB. By 
the computation of stellar wind and atmosphere models for a 0-type star of 60 M© mass, 
Langer et al.l (Il994j ) concluded that the stellar mass lost to the wind bubble is about 32 M© 
in the MS stage, 8 Mq in the Luminous Blue Variable (LBV) stage, and 16 M0 in the WR 
stage. Such a massive wind bubble will produce a significant shell mainly swept up by the MS 
wind at the time of the supern ova or GRB (see FigjID, with a distance from the progenitor 
of about R ~ tens of parsecs (JGarcia-Segura et al.l Il996al : iMirabal et al.l l2003l : iDwarkadas 



20071 ). If the clumps formed during the LBV and WR stages survive lo ng enough, they will 



be hit ting the MS shell roughly at the time of the supernova or GRB (JGarcia-Segura et al. 
1996ah . 



In this paper, we apply the echo emission model to a sample of 36 GRBs with well- 
monitored X-ray light curves by the XRT onboard Swift. We find that most of the light 
curves without significant flares can be reproduced by our model. In this sample, 14 bursts 
have available redshifts by observations. By fitting their X-ray light curves, we find the 
distance R of the wind bubbles from GRB progenitors. A measurement of R might bring us 
closer to finding the mass loss rate, density and wind velocity of the progenitor prior to the 
explosion. 

Our paper is structured as follows: In §21 we first discuss the existence of massive stellar 
wind bubbles around GRBs when they explodes and the subsequent X-ray echo emission. 
We formulate the theoretical light curves of X-ray echo emission in ^ and we introduce 
the temporal break in echo emission and fit the observational light curves of 36 GRBs in ^ 
Discussion and conclusions follow in © 



2. THE ECHO ON THE BUBBLE 



It has been shown by lShao and Pail (120071 ) with a detailed calculation following lMiralda-Escude 



(I1999I ) that, the fiducial flux of the echo emission could be comparable to an already ob- 
served X-ray afterglow at some certain time, given that the prompt X-rays from a GRB are 
significantly scattered by the circumstellar dust. Here, we provide a more concise approach 
by order-of-magnitude estimation. 

Evidence accumulated in the last few years has verified th e long-hypothesized origin of 
soft GRBs in the deaths of WR stars (jWoosley fc Bloomll2006l . and references therein). It is 
generally believed that the progenitor WR star will blow out a hot bubble in the interstellar 
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medium (ISM) by massive stellar wind (JGarcia-Segura et al.lll996al : iDwarkadad 120071 ) . The 
typical radius of the dense shell of swept-up surrounding ISM is of order ten s of parsecs, as 



given by a self-similar solution (jCastor et al.lll975l : iRamirez-Ruiz et al.ll200ll : iMirabal et al. 
2003h 



R 



f 25Mw2 



1/5 



^3/5 



1/5 



100 



M_6.2<5 \ .3/5 



no 



h.6 PC, 



(1) 



where the mass-loss rate during the MS stage M = 10 ^-^ Mq/jy, the typical wind velocity 
V = 10^'^ km/s, and the lifetirne of the MS stage t = 10^'^ yr are adopted for a 35 Mq 



star (JGarcia-Segura et al.lll996bl : iDwarkadad 120071 ) . and the number density of interstellar 
medium n = 1 cm~^ is assumed. 

Correspondingly, the geometry and density profile of a massive wind bubble around a 
WR star are shown in FiglH WR stars - char acterized by high ma s s-loss rate - are believe d 



to form massive dust in their winds (e.g., IWilliams et al.l 119871 : ICrowther 1 12003| . 120071 ). 



Interestingly, all dust shells that were observed were essentially associated with carbon-rich 
WR stars of the latest subtypes of the WC sequ ence only (which are called WC stars herein. 



e.g., lAllen et al.lll972l : iGehrz fc Hackwelllll974j ). WC stars are the final evolutionary phase 
of the most massive (> 4OM0) stars, which, in some sense seem to be favored here as the 
progenitors of some GRBs, since there is where echo emission would be most efficient. 

Assuming that an X-ray burst (i.e., the X-ray counterpart of a prompt GRB produced 
at the final stage of a WC star) with a fiuence Sx is scattered by the dusty shell at a distance 
R with a scattering optical depth Tsca, the echo emission is spread within a characteristic 
time scale td, due to the scattering delay of X-rays with different scattering angles a (see 
FigJUf or the geometry of scatter ing). Since the characteristic scattering angle a is physically 
small (jAlcock fc Hatchettlll978l ). the time scale t^ is consequently much smaller than R/c, 
where c is the speed of light. Therefore, we have a fiducial fiux of the X-ray echo emission 
over a time T^ 



-I echo 






2.5 X 10 



-11 



s- 



X 



2.5 X 10"'^ergscm~2y Vo.l 



n, 



103; 



ergs cm ^ s '^, 



(2) 



where S'x is the total fiuence of the burst observed in X-ray band (say, 0.3-10 keV), the echo 
fiux and the source fiux are assumed to suffer the same interstellar absorption on their way 
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to the observer, and the cosmological fc-correction is ignored. 5'[o.3_io] 
is approximated here, assur aing 5'[n..s-in]/'S'[iR im] = 1/4 and ^[is^iso] 
that 5'[t;n_inn]/5'[ot; t;n] ~ 1 (jSakanioto et aL 



soft X-ray band (jPreece et al 



200C 



= 2.5 X 10"''' ergs cm"^ 

lO^^ergscm"^, given 

20071) and mos t GRBs have a flat spectrum in 



0|; ISakamoto et ahl 120071 ). Note that, a reahstic S'[o.3_io] 



should be essentially larger than the value we adopt h ere, because the X-ray s between [0.3-10] 
keV before and after the BAT trigger interval (Tioo; ISakamoto et al.l 120071 ). especially those 
from the soft tail emission in most GRBs, are not included in this approximation. Since the 
scattering angle a is very small, we have 



td 



;i + z)Ra^ 



10^ 



2c 



R 



100 pc 



T7 S, 



(3) 



where z i s the redshift of the GRB, and a is the typical scattering angle of X-ray echo 
emission JAlcock fc HatchettI Il978l : Ixiose I llQQsh . 



Hence, we have three interesting coincidences: 



1. 



he d usty shell at a distance R of about tens of parsecs assumed by IShao and Dai 
( 20071) is consistent with the massive wind bubb l e around a WR s tar at its final stage 
(JGarcia-Segura et al.lll996al : iMirabal et al.l l2003l : lDwarkadasl 120071 ): 



2. The timescale of the X-ray echo emission on such a wind bubble is about td ~ 10'^ 
s, which i s the typical time s cale when the shal low decay phase dominates in X-ray 
afterglow Jzhang et allbood : iNousek et al.ll2006h : 



3. The X-ray echo flux F ° from such a wind bubble is co nsisten t with most of both 



WQ- S wift and Swift X-ray afterglows around 10^ s (e.g., ICostal Il999l : iNousek et al. 



20061 ). 



These three coincidences motivate the study of echo emission as an additional contribu- 
tor to the X-ray emission. In what follows, we work out the light curve of the echo emission, 
and show the fourth and the most spectacular coincidence that the theoretical light curves 
of X-ray afterglows are consistent with the observed ones. 



THEORETICAL LIGHT CURVE 



Shao and Pail (120071 ) have derived the temporal profile of the echo emission. They 



showed that the feature in the X-ray light curve of a shallow decay followed by a "normal" 
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decay and a further steepening could be reproduced, given that the grains in the dusty 
shell have a distribution over size and that the scattering is depei ident on the w a velen gth of 
X-rays through the Rayleigh-Gans approximation. According to IShao and Pail (120071 ). this 
theoretical light curve can be formulated as 

^echo(^) = A /" ^ /" ^ S{E)T{E,a)j^{x[a{t)]}-dadE, (4) 



S{E)t{E, a)jf{x[a{t)]}-dadE, 

where A is a constant, E is the energy of an X-ray photon, a is the size of a grain in the 
dusty shell, and we define the following functions 

5 



S{E) 



E 



100 keV 



exp 



{5+l)E 

Er, 



(5) 



T{E,a) 



[^ + z)E 
IkeV 



Jiia^j 



x^ 



0.1 /xm 



sm X cos X 



X 



4-5 



(6) 

(7) 



x{a) 



2n{l + z)Eaa 
he 



a{t) 



2ct 



'l + z)R' 



(9) 



wher e 6 and Ep are the spectral index and the observed peak energy of the source (IPreece et al. 




s and Q are t 



le power- law indices in the formula of scattering optical depth (jMauche &: Gorenstein 



Mathis et al. 



Mat his et al. 



19771 : IShao and Pail 120071 ). a_ and a+ are the cutoff sizes of dust grains 



19771 ). [E_, E^] is the energy band of the detector, and h is the Planck constant 



see the Appendix for a detailed derivation) 



4. TEMPORAL BREAK IN X-RAY ECHO EMISSION 



Thanks to the Swift satellite (JGehrels et al.ll2004l ). our interpretation of GRBs has be- 
come more complex. The feature of a shallow decay followe d by a "normal" decay and a fur- 
ther r apid decay of X-ray afterglows appears to be canonical (IZhang et al.ll2006l : iNousek et al. 
20061 ). Coincidentally, we find that this feature can be naturally reproduced by the echo 



emission. 
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Well-monitored X-ray light curves of 36 GRBs (see Fig. [21) were fi tted by our theoretical 
model above (see En. H]), using the online data by lEvans et al.l (120071 ). A total of 14 events 
have known redshift^j (see Tab. [T]) while 22 without known redshifts are fixed to he at z = 1 
(see Ta b. H. We assum e E_ = 0.3 keV and E+ = 10 keV f or Swift XRT JCehrels et al. 



2004: 



Evans et al. 



20071) 



0.005 /im flMathiset al.lll977l ). 6 = 0, and Ep = 200 keV 



(jPreece et al.ll2000l ). The other parameters, s, q, a+, and R are given in Tab. [Hand Tab. 
[21 Apparently, a+ is almost a constant here (i.e., a+ = 0.25 /im, except a^. = 0.5 /zm 



for GRB 050505) and consiste nt with most results from optical observations (JMathis et al. 



19771 : IPredehl fc Schmittlll995l ). 



Here A is a numerical coefficient, which is equal to lA^Ai as shown in the Appendix. 
Al is the coefficient in the integrated GRB spectrum between 0.3 and 10 keV, which is not 
available based on current working frequency of BAT onboard Swift. Similarly, A2 is hardly 
determined by current observations about the X-ray optical depth. A further analysis of the 
parameter Ai and A2 is required once the relevant observational data are sufficient. 

In Fig. [21 all the light curves approach to a steep decay (oc t~^), which is already 
predicted by Eq. (^. Since this is a double integral about E and a, t is an independent 
parameter. Approximately, when x[a(t)] ^ 1, we have 



-lecho 



it) 



OC 



OC 



oc 



FEtit) 

sin x[a{t)] 



x^[a{t)] 
cos'^ x[a{t)] 



cos x[a(t)]^ 1 
x[a{t)] J 1 



OC t 



x^[a{t)]t 
-2 



(10) 



where cos^x[a;(t)] is a high-frequency periodic function of t when x[a;(t)] ^ 1, and will be 
independent on t after the integration over E and a. Therefore, letting a;[a;(t)] = 1, we have 
a critical time 



tc. 



h^cR 



2 



lOS 



E yv ^ 

IkeVy VlOOpc 



1 + z J \0.1/im 
after which, the light curves will be approaching the steep decay (oc t~^ 



;ii) 



^ The redshift data come from http: //www. mpe.mpg.de/$\sim$jcg/grbgen. html 



By fitting the light curves, we can indirectly derive the distant R for a given event at 
redshift z (see Tabled]). These v alues of R appear to be consiste nt with the observed radii of 
wind bubbles around WC stars (JMarston 1119971 : IChu et al.lll999l ). For bursts without known 
redshifts we only list a pesudo value of the bubble radius -Rpseudo assuming a fixed redshift 
z = 1 (see Table [2]). The larger scatter in -Rpseudo values for the latter group appears to be a 
direct consequence of the arbitrary z = 1 prescription. 



5. DISCUSSION AND CONCLUSIONS 



Recently, iLiang et al.l (l2007al ) provided a comprehensive analysis of shallow decays and 
"normal" decay segments in GRBs. They concluded in their sample that most GRBs in 
these two segments have a spectral index (3 around —1.1 (as the flux density Fe oc E^). 
Here we can also get a knowledge of the spectrum of echo emission based on our analysis 
(see Eq. IA9p . Approximately, when x[a(t)] ~ vr at a given time t.,^, we have 



F-^°(t.) ex Flt(t.) 

oc S{E)T{E,a) 
oc E^~\ 



(12) 



where t^, is about an ord er of magni t ude la rger than t^ in Eq. fITT]) . and roughly equals 
most of the t2's defined in lLiang et al.l (J2007al ). Therefore, we should have the spectral index 
/? > (5 — s before tjr, since the spectrum of echo emission is softening in the long run. As shown 
in Fig. [31 the hardness ratio of echo emission is monotono usly decreasin g. Interestingly, it 
reaches a platform (instead of an extremum proposed by iKlose I Il998l) almos t around tjr 
and then keeps decreasing v e ry slowly. If we assume 6 :^ (jPreece et al.ll2000l ) and s ~ 2 
(JMauche fc GorensteinI Il986l : iMeszaros fc Gruzinov 1120001 ), we should h ave >, —2. This 
seems to be inconsistent with the observational result by iLiang et al.l (l2007al ). since they 
ended up with /? ~ —1.1. 

However, we have two reasons that lead us to be content with our result based on current 
analysis. First, since we are dealing with the prompt X -ray emission and the 5 ~ here is the 
result from hard X-rays by BATSE (jPreece et al.ll2000l ). we can imagine that the real average 
S should be larger than 0. Specifically, based on the standard synchrotron emission, 6 should 
equal 1/3, and the synchrotron self-absorption leads to 5 ~ 2. Second, the Rayleigh-Gans 
approximation is adopted here to calculate the differential cross section of scattering. This 
approximation overestimates the echo emission at softer X-ray band (say, in 0.3-1.5 keV), 
where absorption is important, which will change the spectral shape in soft X-ray band and 
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alleviate the dilemma of a dramatically small hardness ratio in Fig. [31 Based on the current 
analysis, it is safe to say that -2 is only a lower limit for (3. Therefore, our result here is 
generally consistent with the observati ons. Thus, t h e echo emission cannot be completely 
ruled out for the reason as discussed by iLiang et al.l (l2007al ) . 



Originally, the opti cal echo emission was used to explain the red bumps observed in 
some optical afterglows (JEsin fc Blandfordll2000l). However , since the typical optical/infrared 
scattering angle would be about a ~ 60° (IWhitd Il979l ). the timescale in Eq.Q would 
be longer, and the optical /infrared echo e mission would be much weaker, and hard to be 
detected ( lReichartll200ll : iHeng et al.l 120071 ). Nevertheless, if the WR wind within a much 
shorter distance [R ~ 10^^ — 10^^ cm) is also taken into account, there might be detectable 
variations in the optical afterglows due to echo emission (JMoran Sz Reichartll2005l ). However, 
the av erage LBV or WR wind bubble a round a massive star at its death would be much 
larger ( Garcia-Segura fc Mac Lowlll995al Jbl). Therefore, the features of echo emission in the 
optical band would be not significant enough, which is possibly the reason that the standard 
external shock model was more successful in the pre- Swift era, when the optical observations 
dominated. 

Furthermore, the dramatic variations in the LBV or/and WR w inds would drive in- 



stabili tie s that produce r adial filaments in the LBV or WR shells (JGarcia-Segura et al. 



1996a 



Qt 



Dwarkadad 



glows (JBurrows et al 



20071) . which could be the origin of some flares in the X-rays after- 



20061 ). Therefore, there might be another timescale due to the angular 



variation of these filaments, which could be given by 



ti 



[1 + z)Ria1 



(13) 



where R{ is the distance of the filaments from the progenitor and ai is the typical angular 
size of the clumps and ripples in the filaments. To account for the mean rat io of the width 
and peak time of the X-ray flares (JBurrows et al.ll2007l : IChincarini et al.ll2007l ). i.e., tf/td ~< 
At/t >~ 0.1, we may have Of ~ 0.1 a ~ 0.1'. 



Sazonov fc Sunyaev I (120031 ) claimed on a different scattering circumstance (molecular 



and atomic matter instead of dust grains) that the opening angle of jets might be an impor- 
tant factor to the echo light curve, based on the popular assumption that the GRB prompt 
emission is beamed. Since only small-angle scattering (a ~ 1') for dust grains is considered 
here, the typical jet properly pointing at us with a much larger opening angle of about few 
degrees will not affect out results. However, the GRB jet might be pointing at us with its 
edge in some case, where our result here would not be valid. Instead, a rapid drop that 
is much steeper than t~^ will emerge at a late time (ISazonov fc Sunvaev II2003I ). This is 



probably the case for some GRBs, e.g., GRB 060526 (JDai et al. 



200 



W 
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in this paper, we have presented the X-ray echo emission in GRB phenomenon and we 
can draw the following conclusions: 

1. The shallow decay phase can still be understood without a dditional contributions from 
the central engine (jZhang et al.l l2006l : iLiang et al.ll2007al ). Interestingly, the feature 
of a shallow decay followed by a "normal" decay and a further rapid decay can be 
reproduced by the echo from the p rompt X-ray emission scattered by the niassive 



wind bubble aro und the progenitor (JGarcia-Segura et al.l Il996al : iMirabal et al.l 12003 



Dwarkadasll2007l ) ■ In general, dust formation seems to be prevailing in carbon-rich WR 



stars of the WC subtype and there is where echo emission would be most efficient. 

2. Echo emission c an be an additional X-ray source that masked the X-ray jet breaks 

( Sato et al.KOOTI') . The r efore, the lack of jet breaks (jPanaitescu et al.ll2006l : lBurrows fc Racusin 



2007 



Panaitescul 120071 : iLiang et al.l l2007bl ) might also be understood via competing 



echo emission. In addition, some results in t he we- Swift era, say, some correlations 



relevant for optical b reaks might still be valid (JFrail et al.ll200ll : iGhirlanda et al.l 12004 



Liang fc Zhandl2005l ) . 



There is an extra temporal break in X-ray echo emission introduced by t^ above. By 
fitting the light curves, we can indirectly measure the locations of the massive wind 
bubbles, which will bring us closer to finding the mass l oss rate, wind velocity, and 
the age of the p rogenitors prior to the GRB explosions (JGarcia-Segura et al.l Il996al : 



Dwarkadasll2007h . 



Even though some features in X-ray afterglows would be better understood if the echo 
emission is taken into account, the external shock model is still the dominant scenario 
for optical/infrared afterglows, since echo emission is negligible in the later. In addition, 
the echo emission model does not completely rule out other interpretations (energy 
injection, etc.) in the framework of the fireball-shock model. Technically, echo emission 
takes place at a distance of ~ tens of parsecs, while the internal or/and external shocks 
are believed to take place within a distance of ~ 10^^ cm. It is more likely that both 
of them are taking place in the X-ray afterglow, which potentially complicates the 
light curves. However, this may raise the problem about GRB radiative efficiencies 
again. Even by assuming that the X-ray emission is the forward shock emission, the 
GRB eff iciency is already t oo high to be accommodated by the leading internal shock 
models (IZhang et al.l 120071 . and references therein). Assuming an echo-dominated X- 
ray emission would exacerbate the situation. 
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A. APPENDIX: DERIVATION OF THE THEORETICAL LIGHT CURVE 



The single-scattering appro ximation is adopted, which is valid up to a scattering optical 



depth, Tsca, of ~0.5 at 1 keV (IPredehl fc Klosd Il996l : iKlose 1 119981 ). Therefore, the echo 



intensity at the Ea rth per unit photon energy per unit size of the grains is (see Fig. 1 of 
Shao and Pail 120071 . for the geometry and quantities defined in this Appendix) 

da 



lt:{OA\t) = FE[t-U{d)]Ta{E)- 



cr=, 



,dfi' 



(Al) 



where FeIP — td(^)] is the flux density of the initial GRB before a time delay t^{9), which is 
also dependent on the viewing angle 6, Ta{E) is the frequency-dependent scattering optical 
depth per unit size of the grains, and the last term is the fraction of the scattered photons, 
which go into the line of sight in a unit solid angle. The zero time point is set as the trigger 
of a GRB on the Earth. is the azimuthal angle along the line of sight. 

Since the initial GRB is usually much shorter than the long-term afterglow, we can 
approximate the initial GRB light curve as a pulse. We also assume that t he initial GRB 
spectrum is given by a Band spectrum in the X-ray band (JBand et al.lll993l ). so we have 



FE[t-t^i9)] = A, 



E 



exp 



{5+l)E' 



x5[t-Ui9)], 



(A2) 



100 keV 

where we only need the segment of the Band spectrum in the lower energ y band. A^ is a 
const ant, and the time delay from a certain viewing angle 6 is given by (JMiralda-Escude 
1999h 



td(^) 



(A3) 



2cDds 

Here the Dirac delta function of t can be translated into a function of 6, using the special 
properties of its own, 

cDds 



s[t-u{e)] 



[1 + z)D.D£ 



6[e-e{t% 



(A4) 
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where the function 9{t) is defined as 9{t) = [2ctDds/(l + ^)-Dd-Ds]^^^ (see the Appendix of 
Shao and DailboOTh . 



We also assume the frequency-dependent scattering optical depth per unit size of the 
grains is 



TaiE) 



A. 



[l + z)E 
IkeV 



0.1 /im 



i-q 



(A5) 



where A2 is a constant with a dimension of \a ^] . This is valid only in a range of the size of 
the grains, say, a G [a_,a+] (jShao and Daill2007l . and the references therein). 



For the last term in Eq. (lAlj) . we have dQ = 6d6d(f), since 6 is very small. We can do 
some rearrangement 

^^ = ^^ ^^sc , (A6) 

(Tscadfi (Tscadr^SC ^Q 

where cXsca is the total scattering cross-section, df^sc = adadcp is defined since a is also 
very small, and we have a = {Ds/Dds)6 in geometry. Neglecting the chemical composition 
or shap e of the grains, the Rayleigh-Gans appro ximation is valid for the differential cross- 
section (jOverbeckl Il965l : lAlcock fc HatchettI Il978l ) 

dc^ _ 2 jf{x) 

where x = 27r(l + z)aEa/hc, and ji{x) = (sinx)/x^ — (cosa;)/x is the first-order spherical 
Bessel function. 

Since the wind bubble is very close to the GRB, we have D^ — D^ and D^s = R- 
Therefore, the echo fiux density at the Earth per unit size of the grains is 



^C(^) 



rccho / 



■), (j)] t)cos6dQ 



ASiE)T{E,a)j!{x[ait)]}^., 



(At 



where A is a numerical constant, and S{E), t{E, a), ji{x), x{a), and a{t) are given by Eqs. 
([5])- ([9]). So, the fiux density and the fiux of echo emission at Earth will be 

Ff^°(t) = r Ff:{t)da (A9) 



and 



F^=^°(t) 



F|^^°(t)dE^ 



(AlO) 



respectively. 
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Table 1. Model parameters of 14 GRBs with known redshifts. 



GRB 


z^ 


Begin (s)^ 


End (s)^ 


s 


q 


a+(/im) 


R{pc) 


xVdof 


GRB 050319 


3.240 


4x 


102 


2 X 10^ 


2.5 


4.4 


0.25 


75 


1.72 


GRB 050505 


4.27 


3x 


103 


2 X 10^ 


2.0 


4.0 


0.5 


50 


1.42 


GRB 050814 


5.3 


1 X 


10^ 


1 X 10^ 


3.0 


4.5 


0.25 


55 


3.53 


GRB 051022 


0.8 


1 X 


10^ 


1 X 10^ 


2.0 


3.1 


0.25 


30 


2.85 


GRB 060210 


3.91 


3x 


10^ 


8x 10^ 


2.0 


3.0 


0.25 


200 


1.78 


GRB 060218 


0.033 


1 X 


10^ 


1 X 10*5 


3.0 


3.1 


0.25 


50 


1.92 


GRB 060502A 


1.51 


4x 


103 


2x 10*5 


3.0 


4.3 


0.25 


65 


1.26 


GRB 060512 


0.4428 


3x 


103 


3x 10^ 


2.0 


3.0 


0.25 


45 


2.63 


GRB 060707 


3.425 


1 X 


103 


2x 10^ 


3.0 


4.5 


0.25 


170 


2.53 


GRB 060714 


2.711 


3x 


102 


1 X 10^ 


2.0 


4.0 


0.25 


95 


2.61 


GRB 060814 


0.84 


8x 


102 


1 X 10^ 


2.0 


4.0 


0.25 


34 


2.33 


GRB 060729 


0.54 


5x 


102 


1 X 10^ 


3.0 


5.5 


0.25 


45 


2.93 


GRB 061110A 


0.758 


6x 


103 


7x 10^ 


3.0 


4.5 


0.25 


76 


2.61 


GRB 061121 


1.314 


2x 


102 


2 X 10^ 


2.0 


4.0 


0.25 


40 


2.30 






















' http : //www . 


mpe . mpg 


de/$\sim$_ 


j cg/grbg 


sn . html 





^Beginning time of the data for fitting 
'^Ending time of the data for fitting 



-17- 



Table 2. Model parameters of 22 GRBs without known redshifts^ 



GRB 


Begin (s)'' 


End (s)'= 


s 


q 


a+(/im) 


-Rpscudo (pc) 


xVdof 


GRB 050712 


4x 


103 


2 


X 10^ 


2.3 


4.3 


0.25 


55 


2.75 


GRB 050713A 


4x 


103 


2 


X 10^ 


2.0 


4.0 


0.25 


45 


2.67 


GRB 050713B 


5x 


102 


1 


X 10^ 


3.0 


4.5 


0.25 


70 


2.21 


GRB 050802 


3x 


102 


1 


X 10^ 


2.0 


3.1 


0.25 


30 


4.64 


GRB 050822 


1 X 


103 


5 


X 10^ 


3.0 


4.5 


0.25 


80 


2.19 


GRB 050915B 


1 X 


103 


5 


X 10^ 


3.0 


4.8 


0.25 


70 


1.60 


GRB 051008 


3x 


103 


4 


X 10^ 


2.5 


3.5 


0.25 


10 


1.93 


GRB 051016A 


4x 


103 


6 


X 10^ 


2.3 


4.0 


0.25 


80 


6.44 


GRB 060204B 


4x 


102 


4 


X 10^ 


2.0 


3.1 


0.25 


25 


1.64 


GRB 060219 


2x 


102 


4 


X 10^ 


2.5 


3.5 


0.25 


45 


1.04 


GRB 060306 


3x 


102 


3 


X 10^ 


2.0 


4.0 


0.25 


50 


1.41 


GRB 060428A 


3x 


102 


3 


X 10^ 


2.7 


4.2 


0.25 


95 


2.26 


GRB 060510A 


1 X 


102 


5 


X 10^ 


3.0 


4.5 


0.25 


10 


2.28 


GRB 060604 


3x 


103 


1 


X 10^ 


2.0 


4.0 


0.25 


80 


1.86 


GRB 060708 


2x 


102 


1 


X 10^ 


2.0 


4.0 


0.25 


50 


2.03 


GRB 060813 


9x 


10^ 


2 


X 10^ 


2.0 


3.1 


0.25 


10 


2.68 


GRB 061004 


3x 


102 


1 


X 10^ 


2.0 


4.0 


0.25 


15 


2.61 


GRB 061222A 


2x 


102 


2 


X 10^ 


2.1 


4.2 


0.25 


30 


1.90 


GRB 070129 


1 X 


103 


2 


X 10^ 


3.0 


4.5 


0.25 


90 


1.63 


GRB 070328 


1 X 


102 


6 


X 10^ 


1.9 


4.0 


0.25 


5 


4.75 


GRB 070419B 


4x 


103 


4 


X 10^ 


2.0 


3.1 


0.25 


20 


2.74 


GRB 070420 


3x 


102 


4 


X 10^ 


3.0 


3.1 


0.25 


15 


3.07 



'^We assume z=l to carry out the fitting 
^Beginning time of the data for fitting 
•^Ending time of the data for fitting 
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Fig. 1. — Geometry and density profile of a massive wind bubble around a WR star. The 
major scattering takes place at the massive and dense shell (marked with brick- wall pattern), 
which is a main dusty shell around the WR star (or more probably a WC star). Inside the 
shell, the interior structure of the bubble (filled with light gray) is ignored, which is mainly 
determined by the LBV and WR winds and not important for scattering, since the density is 
much lower there and essentially no dust survives close to the WR star. Outside the shell, the 
interstellar medium (ISM) dominates (filled with heavy gray). Meanwh ile, the gas density p 



is also plotted with logar ithmic scale in dot-dashed line (not in scale, see lGarcia-Segura et al. 



1996a 



Dwarkadad 120071 ). 
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Fig. 2. — The X-ray light curves of 36 GRBs observed by Swift XRT. WT mode is in gray, 
PC mode is in black. The solid lines are the theoretical light curves given by our model (Eq. 
H]). We have E_ = 0.3 keV and E^ = 10 keV for Swift XRT, and we assume a_ = 0.005 /im, 
6 = 0, Ep = 200 keV. The other parameters, s, q, a+, and R are given in Tab. [Hand Tab. O 
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Fig. 3. — Evolution of the hardness ratio of echo emission with different shell distance R. 
We assume z = 1, q = 3.5 and a+ = 0.025 //m. Other parameters are the same with Fig. 
[2j The dramatically small ratio could be alleviated only if the absorption in 0.3-1.5 keV is 
considered. 



